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■ Abstract: We consider a generalization of Einstein's general theory of relativity such 
psj ' that it respects local scale invariance. This requires the introduction of a scalar and a 

■qj ■ vector field in the action. We show that the theory naturally displays both dark energy 

I and dark matter. We solve the resulting equations of motion assuming an FRW metric. 

! The solutions are found to be almost identical to those corresponding to the standard 

! ACDM model. 

C^: 1 Introduction 

r-| ' In recent papers [1, 2] we have studied a scale invariant extension of the general theory of 
relativity. The theory postulates a scalar field $ and a modification of the gravitational 
action such that [3,4], 

■ where G is the Newton's gravitational constant, R is the Ricci scalar and /3 is a dimensionless 
constant. The modified action involves no mass parameter and is invariant under the scale 

QQ I transformation. It is also convenient to introduce the concept of pseudo-scale invariance. In 
four dimensions the pseudo-scale transformation can be written as [3,4], 



X X . 



O 

> 

X 

$ ^ -K/A, 

^ ^ ^/A3/2. (2) 

where x is the space-time coordinate, g'^" the metric, a vector field and \I' a spin half 
field. The scale transformation can be expressed as a combination of the pseudo-scale and 
the general coordinate transformations. Hence as long as general coordinate invariance is 
respected, pseudo-scale invariance is equivalent to scale invariance. 
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We also considered a theory with local scale or pseudo-scale invariance [3-12]. The 
transformations in this case are given by Eq. with the parameter A — > A(x). In this case 
we need to introduce the Weyl vector meson, S^, besides the scalar field. Under pseudo-scale 
transformation the vector field transforms as [3,4] 

S^^S^- jd.lniM^)) ■ (3) 

The Lagrangian in this case, ignoring all other fields besides the $ and S*^, may be written 
as 

C = -^<^^R + Cmatter (4) 
o 

where 

Cmatter = ^9^' {V ,<i>) {V,<^) - - ^g'^" g'^'^ E , (5) 

/ is the gauge couphng constant, E^^ = d^S^ — dyS^, R is the modified curvature scalar [3,4], 
invariant under local pseudo-scale transformation, and 

^ d,-fS, (6) 

is the gauge covariant derivative. The scalar R is found to be 

R = R- 6fS% - 6/^5'^^'^ (7) 

where R is the standard curvature scalar. The mass of the vector boson 5*^ has been con- 
strained by cosmological observations to be very light, less than 400 eV, or very heavy, 
greater than Planck Mass [13]. 

A theory with local scale invariance is aesthetically pleasing since it does not contain 
any mass parameter in the action. This also implies that we cannot add the cosmological 
constant term in the action. Hence as long as scale invariance is unbroken, cosmological 
constant is identically zero. In Ref. [1,2] we argued that scale invariance is broken by the 
cosmological time evolution or equivalently by initial conditions. This phenomenon is called 
cosmological symmetry breaking in Ref. [1,2]. The authors argue that the universe is a time 
dependent solution of the equations of motion. All phenomena are described by making a 
quantum expansion around this time dependent solution. Hence if the symmetries of the 
action are not respected by this solution then these symmetries will be hidden, in analogy 
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to the phenomenon of spontaneous symmetry breaking. However as shown in Ref. [1,2], 
cosmological symmetry breaking is not the same as spontaneous symmetry breaking. 

Once scale invariance is broken cosmologically, the theory generates the dimensionful 
parameters such as the Newton's constant and the vacuum (or dark) energy [14-23]. These 
essentially get related to the initial conditions imposed on the scalar field. The fact that this 
theory leads to dark energy has also been noticed in Ref. [24]. However the precise nature of 
identification is different from that in Ref. [1,2]. In Ref. [1,2] the authors speculated that if 
scale invariance is broken cosmologically, then this symmetry may not be anomalous. This 
symmetry may also tame the quantum corrections to the vacuum energy, hence avoiding 
the fine-tuning problems in cosmological constant [25-31]. Alternate proposals to solve the 
cosmological constant problem are discussed in Ref. [25,32-39]. 

In the case of global scale invariance, the authors in Ref. [2] found a solution with $ equal 
to a constant. In this case our theory reduces to a scalar-tensor model. The cosmological 
implications of such models have been studied extensively in recent literature [40-42]. The 
value of the constant field $ can be chosen to fit the experimental value of the gravitational 
constant. For local scale invariance, $ is constant in the gauge 5*0 = 0. In Ref. [2], the 
authors set Si = 0, where the index i = 1,2,3. However in general Si is not zero. In the 
present paper we give a general solution to the equations of motion of this model. We find 
that the solutions naturally lead to both dark energy and dark matter. 

In Ref. [1,2], the authors had suggested that the scalar field $ may be the Higgs mul- 
tiplet. They argued that cosmological symmetry breaking also leads to a breakdown of the 
electroweak symmetry. In this theory the Higgs particle is absent from the particle spectrum 
and acts as the longitudinal mode of the vector field S*^ [3,4]. Although this is a consis- 
tent picture, the field $ may also be associated with a scalar field in a grand unified theory 
(GUT). In the present paper, we assume $ to be a real scalar field. The generalization to the 
case where $ may be a standard model or a GUT scalar field multiplet is straightforward. 

2 Equations of motion 

The Einstein's equations and the equations of motion for $ and S^ following from Eq. (jlj) 
are given in Ref. [2]. Here we assume that, at leading order, all the fields are independent of 
space coordinate and depend only on time. In this case the equations simplify considerably. 
We display these equations below, correcting some typographical errors in Ref. [2]. The 
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time-time and space-space components of the Einstein's equation are 



(8) 



and 



(9) 



respectively. Here dots represent derivatives with respect to time. The equation of motion 
for the scalar field is 



• a 



a P 



$ - f<^So + 3$- - fS.S^'<i> + A<l>^ - 3/<l>5'o- + ^<^R = 



a 4 



The corresponding equations for 5*0 and Si are 



(10) 



a 2 







(12) 



respectively. Since the theory has local scale invariance, we need to fix the gauge in order 
to obtain a unique solution. We choose the gauge 5*0 = 0. In this case Eq. ffTTj) imphes that 
the scalar field $ is a constant, independent of time. 

Thus, we set $(t) = t] = constant with the above mentioned gauge choice. The nonzero 
constant value of t] essentially acts as dark energy in the present model. The existence of 
dark energy is naturally implied by cosmologically broken scale invariance. The resulting 
0-0 and i-j components of the Einstein's equations can then be written as. 



/3 



X 



Sf 

:^^+2^ 



+ 1 + 



3/3 \ fWSf 



2a? 



(13) 



and 



a a 



37]' 2- + 



a 



4 



3X 



^ / 3^\ Pv'Sr 
2a^ \ 2 2a2 



(14) 
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respectively. In Eqs. [12] and [HI sum over the subscript i is implied in terms containing Sf 
or S^. The equations for i] and Si become, 

4 . ^^^^ 





l + f)/^^^5. = (16) 
respectively. 

The background constant value of the field $ can be related to the Planck mass, Mp, by 
the relation [1-4], 

Pv^ - ^ (17) 

The mass of the vector field, S^, is found to be, 

M|=(l + f)/V. (18) 
If the vector field 5*^ = 0, then the Hubble constant, Hq, is given by, 

[Y IT Mp 

Finally the vacuum energy density is given by, 

PA = (20) 

From the relationship between the Hubble constant and the Planck mass, it is clear that 
either A is extremely small or (3 is extremely large. If we assume (3 ~ 0(1), then A is 
found to be of order 10"^''. In Ref. [1,2] the authors argued that this small value by itself 
may not lead to fine tuning problems. These would arise if the quantum corrections at 
each order require very delicate cancellations to maintain the small value of this parameter. 
The quantum corrections in this model have so far not been computed. It is of course also 
important to explain the origin of such a small parameter. This issue is beyond the scope 
of the present paper. However we speculate that such small values of scalar coupling might 
arise due to the well known triviality of scalar field theories [43,44]. The continuum scalar 
field self coupling is driven to zero by renormalization group analysis. However at length 
scales smaller than the Planck length, it may not be appropriate to treat space-time as a 
continuum. The discrete nature of space time might generate a small value for the scalar 
self coupling. 
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3 Cosmological Solution 



The equation of motion for Si, Eq. (1161) . is similar to that of a damped harmonic oscillator, 
with weakly time-dependent frequency and decay terms. We seek a solution of the form 

S^ = niS (21) 
where is a constant unit vector. The solution for S can be expressed as, 

S = i?e(^e"^^^^~*^^i^^ + 5e"^^'^* + *^^i'^^l , (22) 



where, A and B are assumed to be slowly varying functions of time and = u"^ — ^ , u"^ = 
(l + ^) Z^??^ = Mg and H = a/a. By substituting this solution in Eq. (fT6l) and neglecting 
second derivatives of A and B, we get. 



(23) 



where, fci and k2 are constants of integration and are, in general, complex. Since uj » H 
we see that A and B vary very slowly with time compared to other terms in S^. The most 
rapidly varying terms are those containing / ujidt in the exponent. Due to these terms. Si 
fluctuates rapidly with time. 

Due to the presence of a vector field Si in our theory, our cosmological solution naturally 
contains a constant three dimensional unit vector defined in Eq. |2U This vector defines a 
direction in space and hence breaks rotational invariance. However the background metric is 
still isotropic since the vector does not contribute to Einstein's equations. Furthermore it 
is unlikely to lead to very large observable consequences of the breakdown of isotropy. This 
is because the field S^ does not directly interact with visible matter [3,4]. Nevertheless, it is 
extremely interesting to determine the cosmological predictions of this breakdown of isotropy 
in view of several observations which indicate a preferred direction in the universe [45-52]. 
Models in which vector fields acquire nonzero vacuum or background values have also been 
considered by many authors [53-68]. It has been argued that many of these models, which 
lead to prolonged anisotropic accelerated expansion, are unstable [69]. In our model the 
vector field does not directly lead to anisotropic expansion, even though it acquires a non- 
zero background value. 
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3.1 Leading Order Solution 

At leading order we can assume that A and B are time independent. The leading order 
solution for S can then be written as 

S^^{A'cos9 + B'sm9) (24) 
V o 

where 9 — J uidt and A' & B' are some real constants. We also find 

1 / H 



where, p = A' cos 9 + B' sin 9, q = —A' sin 9 + B' cos 9. We next substitute these into the 0-0 
component of the Einstein's equation. We define, 

^^.^t^^ul^Sl (26) 

This essentially acts as the contribution to the energy density provided by the field S^. We 
find, 

1 r _ _ _ F7-2 _ _ 1 

(27) 



1 

PSi 



2a3 

where, p^ + q^ = A'"^ + 5'^, which is a constant. Since uji is very large. Si is a rapidly 
oscillating function of time. Hence it is reasonable to replace the oscillatory functions with 
their time averages. After averaging over time, (p^ — g^) ^ and pq 0. Hence, 

P^. = ^ + (28) 
We next consider the i-j component of the Einstein's equation. We define, 

which effectively acts as the contribution of the Si field to pressure. After substituting the 
time averaged values for the oscillatory functions, we get 



Si 



(30) 



Hence, we find that the field Si essentially acts as the cold dark matter. Its energy density 
ps^ varies as and its pressure Ps^ is zero at leading order. A similar phenomenon is seen 
in the case of coherent axion oscillations [70-74] . 
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The modified Einstein's equations, at leading order, can now be written as 

- If"^ - 3i^(^^ - (31) 

•■ • 2 \ 

CL CL A r, . , 

2- + - = • 32 



Eq. (13111 generalizes the expression for the Hubble constant, Eq. (jT9l) . for the case when the 
vector field is non-zero. 



3.2 Corrections to the leading order 

We next calculate the corrections to the leading order result by taking into account the time 
dependence of the coefficients A and B in Eq. fl22|) . Substituting for A and B, from Eq. 
in Eq. we find, 



where, 



S 
S 



1 



^yuJla 



[Q cos{6 - x) + P sm{6 - x)] 



1 



y/iUia 



H f X 



2 



l]U +{uJi-x)V 



U = N cose + M sine ,V = -N sine + M cos9 
and M = P cos x + Qsinx ,N = —P sinx + Q cos x . 



(33) 
(34) 



Here, x 



sm 



-1 H _ I 



cos — , X = H/AuJi = —LOi/H and P Q are some real constants. 



Substituting these in Eq. 

1 



we get. 



(U' + V')iu' + — (U' -V')~iUiH [1 UV 

4 \ uJiJ 



X y 



+ (x^ — 2uJix) 



'IP_ 

M 



(35) 



The third term on the right hand side reduces to uiHUV, since x/ui << 1. The fourth 
term simplifies to —HV^/2 , if we neglect terms suppressed by factors of H /2uji. Hence we 
find 



Ps, 



2ujia^ 



4 



(36) 
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We again substitute time averaged values for rapidly oscillating functions. This sets ([/^ — 
V^) -^0,UV^O and {U^ + V^) = (M^ + A^^) = + Q\ which is a constant. A leading 
order expression for H can be computed using Eq. ( pTl) . We find 



H = --^JA'^ + B'^)uj\ (37) 



Thus, we get, 



2a3 48/3a3M5 6/5r72a6 ' ^ ^ 

The leading term varies as as already found in the previous section. Here, we also find 
two subleading terms. One of these falls as and the second falls much faster, as a ^, as 
the universe expands. We similarly find the corrections to the pressure term Ps-. We find 
that, using Eq. ((291), 



3P. ^ 



{U'-V'Pl + u.HUV + ^V^ 



Again, substituting time averaged values for the rapidly oscillating functions, we get 



(39) 



v+o')=^^^i^i^;i^. (40) 



Hence, we get a small correction term to P^., which also decays rapidly as as the universe 
expands. 

The 0-0 and i-j component of the Einstein's equations can, now, be written as, 

gg„^H^ = V+'^' + f P-^) (41) 



and 



4^ l a + a^i 4^ + 2a;ia3 4 ^ ^ 



respectively. The first of the above two equations can be cast in the form, 

1 = fiA + ^s. (43) 

where 

O - O - 

"A — ; "5i — , 

Per Per 
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PA = -V ,PS.= ^^^^3 -j) and p.. = -r^ i7 . 

Eq. (ia]) looks like the ACDM model with = ^s,- 

Thus, the energy density p^. and the corresponding pressure P^. of the vector field Si 
including the correction terms, are obtained as 

Cl / 2 C2 
Ps, = 7^ J + — 



and 



C1C2 



6 



where. 

In the limit x ^ and tui — tu, we find A' = Q/ ^/uj and P' = P/ y/uj. 

We can make an estimate of the term (P^ + Q"^). Since the recent cosmological observa- 
tions support a fiat ACDM model, we can equate the second term on right hand side of Eq. 
( HTj) . evaluated at present time, to pufi- The contribution due to H is negligible. Thus, we 
get, 

where Qm is computed at the current time. 



4 Including the contribution due to radiation 

In this section, we obtain a set of dynamical equations to study the evolution of different 
components of the universe since the beginning of the radiation dominated era. For this pur- 
pose we introduce, by hand, the contribution due to radiation. We expect to reproduce the 
usual Big Bang evolution where radiation dominates at early times, followed by dark matter 
and dark energy dominated eras, respectively, at late times. We introduce a radiation term 
with energy density pn and it's corresponding pressure term Pr in the energy- momentum 
tensor T^^. The resulting equations are solved numerically. 
It is convenient to introduce the following variables. 
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3/3 a'^T]^ ^ ' 
2 / 3/?\ PS^ _ 

3/3 a ~ 

Here, fig- = f2i + f22, p_R,o is the radiation energy density in the current era and the prime 
denotes derivative with respect to In a. Hence for any function /, 

df Idf 
din a H dt 

With these variables, we can cast the equations (fTSll . (fTSl) and ( fT6l) . along with pn and Pr, 
in a dimensionless form, to obtain the following set of equations : 

X' = X{2 - _ ^2^ ^ ^^^^ 

Y' = -Y{2X^ + Z^) - '^XZ , 

Z' = Z(1-2X^ - Z^) + -XY , 



2 



R' = -2R{2X^ + 



where, k = ^^f-^ = VSMpMs/ y/2¥^. 

We studied the dynamical equations numerically from the beginning of radiation dom- 
inated era (In a = —29) till today (In a = 0). The results are presented in Fig. ([T]). In 
the graphs we only show results for the range Ina = [—14, 0], as radiation is the only dom- 
inant component in the omitted regions. The plots show the results for three values of 
K = 50, 200, 500. The initial conditions for these three cases have been chosen so as to match 
the final observed values of Qm and Qa [75-78]. 

As is evident from the plots, varying k, varies the frequency of oscillations. Besides that 
the results are almost identical, as long as k >> 1. This can be understood from the 
expression of k. For fixed values of Mp and pv, increasing k, increases u or M5, which 
implies more rapid oscillations. Furthermore as seen from our analytic results, applicable 
when radiation energy density is negligible, we reproduce the standard ACDM model in the 
large k limit. 
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Figure 1: The ratio of energy density to the critical energy density, Qi, for different compo- 
nents function of ln(a) for k 



50,200,500. 
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5 Conclusions 

We have analyzed a locally scale invariant generalization of Einstein's gravity. The theory 
requires introduction of a scalar and a vector field. The scale invariance in the theory is 
broken by a recently introduced mechanism called the cosmological symmetry breaking. We 
have shown that this theory naturally leads to both dark energy and dark matter. Due to 
scale invariance the cosmological constant term is absent in the action. The solutions to 
the equations of motion admit a constant, non-zero value of the scalar field, which leads 
to a small cosmological constant or dark energy. The cold dark matter arises in the form 
of vacuum oscillations of the vector field. We have shown that the theory behaves very 
similar to the ACDM model with negligible corrections. The precise values of the energy 
densities of different components are fixed by the initial conditions. Some of the parameters 
in the model take very small values and it is necessary to find an explanation for such small 
values. Furthermore it is important to compute quantum corrections in this model since 
that will determine whether the model suffers from fine tuning problems. The model can 
be generalized to include the standard model fields. The scalar field may then be identified 
with the Higgs multiplet. In this case the Higgs particle is absent from the particle spectrum 
and hence provides a very interesting test of the model. Alternatively the scalar field might 
be identified with a GUT scalar field multiplet. This possibility has so far not been studied 
in the literature. 
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